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In the present paper we study the existence, uniqueness and asymptotic stability of 
the T-periodic solutions for the system 

x = eg(t,x,e), (1) 

where e > is a small parameter and the function g £ C°(R xl'x [0, 1], R k ) is T-periodic 
in the first variable and locally Lipschitz with respect to the second one. As usual a key 
role will be played by the averaging function 

T 

g {v) = J g{r,v,0)dr, (2) 
o 

and we shall look for those periodic solutions that starts near some v £ ^ 1 (0). 

In the case that g is of class C 1 , we remind the periodic case of the second Bogolyubov's 
theorem ([2], Ch. 1, § 5, Theorem II) which represents a part of the averaging principle: 
det (go)'(vo) assures the existence and uniqueness, for e > small, of a T-periodic 
solution of system (J\) in a neighborhood ofvo, while the fact that all the eigenvalues of the 
matrix (goY^o) have negative real part, provides also its asymptotic stability. 

It was Mitropol'skii who noticed that various applications require the generalization of 
the second Bogolubov's theorem for Lipschitz right hand parts. Assuming that g is Lip- 
schitz, g £ C 3 (M fc ,lR fe ) and that all the eigenvalues of the matrix (^o)'^) have negative 



real part Mitropol'skii developed an analog of the second Bogolyubov's Theorem proving 
the existence and uniqueness of a T-periodic solution of system (DO) in a neighborhood 
of Vo- There was a great progress weakening the assumptions of Mitropol'skii in his ex- 
istence result (see Samoylenko [H] and Mawhin [8|) and some analogs of his uniqueness 
result for equations with monotone nonlinearities have been proposed (see papers by A.I. 
Perov, Yu.V. Trubnikov, V.L. Hackevich [10], [32]). Nevertheless, the asymptotic stabil- 
ity conclusion of the second Bogolyubov's Theorem remained to be not generalized for 
Mitropol'skii's settings (namely, when g is Lipschitz) for a long time. It has been done 
recently by Buica-Daniilidis in [3] for a class of functions v i— ► g(t,v,0) differentiable at 
Vq for almost any t G [0, T], but it is assumed in [3] that the eigenvectors of the matrix 
(9o)' ( v o) are orthogonal and that the function g has a continuous Clarke differential, that 
is not easy to check in applications. 

In the next section of the paper assuming that g is piecewise differentiable in the second 
variable we show in Theorem 2 that Mitropol'skii's conditions imply not only uniqueness, 
but also asymptotic stability of a T-periodic solution of system (pQ) in a neighborhood of 
vo. In other words we show that Bogolyubov's theorem formulated above is valid when 
g is not necessary C 1 . Theorem 2 follows from our even more general Theorem 1 whose 
hypotheses do not use any differentiability neither of g nor of g . 

1. Throughout the paper Q C M. k is some open set and for any 5 > the 5-neighborhood 
of v G M. k is denoted by Bg(v ) = {v G M. k : \\v — v \\ < 5}. We have the following main 
result. 

Theorem 1. Let g G C°(R x Q x [0, and Vq G f2. A ssume the following four 

conditions. 

(i) For some L > we have that \\g(t, Vi, e) — g(t, v 2 , e)\\ < L\\vi — v 2 \\ for any t G 



[0,T], Vl ,v 2 eQ,ee [0,1]. 



(it) For any 7 > there exists 5 > such that 



J g(r,v 1 + u(T),s)dr - J Q g(T,v 2 + u(r),s)dr 

-Jo J^i,0)(iT + J T g(r,v 2j 0)dr < ^\\vi - v 2 \\ 



for any u G C°([0, T], R k ), \\u\\ < 5, v u v 2 G B s {v ) and e G [0,5]. 

(Hi) Let g$ be the averaging function given by (d[j and consider that go(vo) = 0. 

(iv) There exist q G [0,1), a, 5o > and a norm || • ||o on M. k such that \\vi + ago(vi ) 
-v 2 - ag (v 2 )\\ < q\\vx - v 2 \\ for any v u v 2 G B So {v ). 

Then there exists 5\ > such that for every e G (0, 5\] system (J\) has exactly one T- 
periodic solution x £ with x e (0) G Bs 1 (vq). Moreover the solution x £ is asymptotically stable 
and x £ (0) —> v as e — ► 0. 

When solution x(-, v, e) of system ([Tl) with initial condition x(0, v,e) = v is well defined 
on [0,T] for any v G Bg (v ), the map v x(T,v,e) is well defined and it is said to be 
the Poincare map of system ((]]). The proof of existence, uniqueness and stability of the 
T-periodic solutions of system (OQ) in Theorem 1 reduces to the study of corresponding 
properties of the fixed points of this map. More precisely, in order to prove Theorem 1 we 
represent x(T, v, e) as 



Then we show that the function g £ satisfies a Lipschitz condition uniformly in e > in 
the ball Bs (v ), where 5 > is a fixed small constant and moreover that g £ satisfies the 
following analog of the property (ii): for any 7 > there exists 5 G [0, 5q] such that 



T 








9e(Vl) ~ 9o(Vi) -gefa) +£002)11 < 7lK ~V 2 



for all vi, v 2 G Bg (vq) and e G [0, 5}. 

This allows to conclude that if / + eg does not exceed 1 — eq (where q is a constant) 
then the Lipschitz constant of the function I + eg £ does not exceed 1 for e > sufficiently 
small. Thus the conclusion of theorem 1 follows by applying the asymptotic stability of the 
Poincare map fixed point lemma (see e.g. [6], lemma 9.2). It turns out that the Lipschitz 
constant of the function I + ego does not exceed 1 — eq indeed, where q = (1 — q)/a, 
providing that the assumptions of theorem 1 are satisfied. 

In general it is not easy to check assumptions (Ej) and (Ev]) in the applications of 
Theorem 1. Thus we give also the following theorem based on Theorem 1 which assumes 
certain type of piecewise differentiability instead of (Jn|) and deals with properties of the 
matrix (go)'(^o) instead of the Lipschitz constant of go. For any set M C [0,T] measurable 
in the sense of Lebesgue we denote by mes(M) the Lebesgue measure of M. 

Theorem 2. Let g G C°(M x Q x [0,1], R k ) satisfying Let g be the averaging 
function given by $M) and consider Vq G Q such that go(vo) = 0. Assume that 

(v) given any 7 > there exist 5 > and M C [0, T] measurable in the sense of Lebesgue 
with mes(M) < 7 such that for every v G Bj(vq), t G [0, T] \ M and e G [0,5] we 
have that g(t, -,e) is differentiate at v and \\g' v (t, v,e) — g' v (t,Vo, 0)|| < 7. 

Finally assume that 

(vi) go is continuously differentiable in a neighborhood of Vq and the real parts of all the 
eigenvalues of (go)'(vo) are negative. 

Then there exists 5\ > such that for every e G (0, Si], system (CP has exactly one T- 
periodic solution x e with x £ (0) G B^iyo). Moreover the solution x £ is asymptotically stable 
and x £ (0) —> v as e — ► 0. 



For proving Theorem 2 we observe that the property (v) implies (ii), while the prop- 
erty (vi) implies (iv). The latter observation is totally based on the lemma 2 from the 
Krasnosekskii's book [5] (P. 91). 

2. Consider now an application of the result announced. 

In his paper [I] Hogan first proved the existence of a limit cycle for the nonsmooth van 
der Pol equation + — l)u + u = 0. By means of the results from the previous section 
we derive now the resonance curves determining the location of stable periodic solutions 
of the perturbed equation 

u + e (\u\ — 1) ii + (1 + ae)u = e\ sint, (3) 

where a is a detuning parameter and eAsini is an external force. 

Following Andronov and Witt [1] we are concerned with the dependence of the ampli- 
tude of 27r-periodic solutions of (3) on a and A providing that e > is sufficiently small. 
By the other words we look for such A and (ft which determine 27r-periodic solutions of (3) 
converging as e — > to 

u(t) = A sin(t + 0), u(t) = A cos(t + 0). (4) 

The assumption (v) of theorem 2 applied to (3) (after rewritting (3) in the 
form (1)) turns out to be satisfied for any a, A G E, while the assumption 
det(go)'((Asm<f), Acoscfi)) ^ (required for validity of (vi)) brings us to the following 
algebraic equation 

A, (» , + ( 1 -^i)') =A ' (5) 

Furthermore the negativity assumption for the real parts of eigenvalues of (go)'(M, N) 



takes the form of the following two inequalities 

39 

vr 2 (l + a 2 ) + — (M 2 + N 2 ) - 4ttVM 2 + N 2 > 0, (6) 
9 

and 

2 (tt - 2v / M 2 ~TiV 2 ) < 0. (7) 

Therefore theorem 2 allows to conclude that the solution u{t) = A sin(t + 0) generates 
asymptotically stable 27r-periodic solutions of equation (3) providing that A, M = A sin 0, 
N = Acoscj) satisfy (5)~(7). 

We note that the corresponding assumptions for the classical Van der Pol equation 
u + e {u 2 — 1) u + (1 + ae)u = eX sin t have the following form (see the formulas (5.21) and 
(16.6)^(16.7) of Malkin's book [7]) 

a, ("' + ( 1 -t) , )= a '- 

1 + a 2 - (M 2 + iV 2 ) + — (M 2 + N 2 ) 2 > 0, 

16 

2 — (M 2 + iV 2 ) < 

respectively. 
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